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Interstellar scintillations of pulsars PSR B0809+74 and B0950+08 have been studied using obser- 
vations at low frequencies (41, 62, 89, and 112 MHz). Characteristic temporal and frequency scales 
of diffractive scintillations at these frequencies have been determined. The comprehensive analysis 
of the frequency and temporal structure functions reduced to the same frequency has shown that 
the spectrum of interstellar plasma inhomogeneities toward both pulsars is described by a power 
law. The exponent of the spectrum of fluctuations of interstellar plasma inhomogeneities toward 
PSR B0950+08 (n = 3.00 ± 0.05) appreciably differs from the Kolmogorov exponent. Toward 
PSR B0809+74 the spectrum is a power law with an exponent n = 3.7 ± 0.1. A strong angular 
refraction has been detected toward PSR B0950+08. The distribution of inhomogeneities along the 
line of sight has been analyzed; it has been shown that the scintillations of PSR B0950+08 take 
place on a turbulent layer with enhanced electron density, which is localized at approximately 10 pc 
from the observer. For PSR B0809+74 the distribution of inhomogeneities is quasi-uniform. Mean- 
square fluctuations of electron density on inhomogeneities with a characteristic scale po ~ 10^ m 
toward four pulsars have been estimated. On this scale the local turbulence level in the 10-pc layer 
is 20 times higher than in an extended region responsible for the scintillations of PSR B0809-I-74. 



1. INTRODUCTION 

Pulsars represent a good tool for the study of in- 
terstellar plasma (ISP), because they possess very 
small angular sizes and intense emission. The in- 
vestigation of their intensity fluctuations in the 
frequency-time domain allows us to study the ISP 
spectrum in various directions of our Galaxy. As 
was shown in [1, 2], the observational data on pul- 
sar scintillations are statistically well described by a 
power-law spectrum of inhomogeneities with an ex- 
ponent n — 3.67 (i.e., by the Kolmogorov spectrum) 
in a very broad range of spatial scales p. However, a 
scatter in the values of the spectrum exponent n is 
observed for the same value of p for different sources; 
this testifies that in local directions of the Galaxy the 
spectrum can differ from the Kolmogorov spectrum. 
The application of a new method for the study of ISP 
using a comprehensive analysis of a structure func- 
tion obtained from multifrequency observations of 
pulsar scintillations has revealed a difference of the 
spectrum from Kolmogorov's for PSR B0329-t-54, 
B0437-47, and B1642-03 [3-5]. The purpose of 
this paper was a study of the ISP spectrum ISP 
toward nearby and powerful at meter wavelengths 
pulsars PSR B0809-h74 and B0950-)-08. The dis- 
tances to these pulsars and their velocities are known 
from parallax measurements [6]: R — 433 pc, V — 
102 km/s and R = 262 pc, V = 36.6 km/s for 
PSR B0809-t-74 and B0950-h08 respectively. Since 
they are the nearest pulsars, the analysis of their 
scintillations allows us to study the region of inter- 
stellar plasma nearest to the Sun. 

The obtained earlier observational data demon- 
strate the presence of three components of interstel- 



lar turbulent plasma in the solar neighborhood at 
distances of the order of 1 kpc and less. The first 
component is turbulent plasma with a statistically 
quasi-uniform distribution in the space between spi- 
ral arms (component A in the classification of [7, 
8]). The second component represents a cavern with 
a depleted electron density (rte — 0.005 cm"'^) in the 
solar neighborhood with a scale of about 200-300 pc 
perpendicularly to the Galactic plane and 50-100 pc 
in the Galactic plane. This cavern has been detected 
in X-ray observations [9] . The presence of a cavern 
with a lowered turbulence level on a scale of the or- 
der of 100 pc has been revealed also by the analysis 
of scintillations of pulsars with dispersion measures 
from 3 to 35 pc cm~'^ [10]. The third component 
corresponds to a layer with an enhanced turbulence 
level at a distance of about 10 pc from the Sun. 
This layer has been detected in observations of in- 
terstellar scintillations of quasars in the centimeter 
waveband [11, 12]. In [5] this layer was classified as 
component C. In [5] it was also shown that the pul- 
sar PSR J0437— 4715 scintillates on inhomogeneities 
of layer C. 



2. OBSERVATIONS AND DATA 
PREPROCESSING 

The observations of PSR B0809-h74 and 
B0950-I-08 were carried out on the BSA and 
DKR-1000 radio telescopes of the Pushchino Radio 
Astronomy Observatory (Astro Space Center, 
Lebedev Institute of Physics, Russian Academy 
of Sciences) at frequencies 41, 62.43, 88.57 MHz 
(DKR-1000), and 111.87 MHz (BSA) in December 
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Parameters 



of 



the 



observations 



Parameter 


/, MHz 


41 


62.43 


88.57 


111.87 




PSR B0809+74 




At, ms 


5.12; 2.56 


5.15 


5.38 


2.56 




128; 64 


96 


96 


128 


A J" 1 TT 

A/, kHz 


1.25; 20 


20 


20 


20 


T, s 


51.68 


51.68 


38.76 


19.38 




PSR B0950 + 08 




At, ms 


5.12; 2.56 


2.56 


2.56 


2.56 


iVeh 


64 


96 


128 


128 


A/, kHz 


1.25; 20 


20 


20 


20 


r,s 


20.24 


15.18 


20.24 


5.06 



2001 - January 2004. As the interference situation 
at low frequencies was complex, we have used for 
the analysis only those records where interference 
was small. We received linearly polarized emis- 
sion. Two multichannel receivers were used: a 
128-channel receiver with a channel bandwidth 
A/ = 20 kHz at frequencies 88.57, 62.43, and 
111.87 MHz as well as a 128-channel receiver with 
a channel bandwidth A/ = 1.25 kHz at a frequency 
of 41 MHz. The time of observations on BSA 
(frequency 111.87 MHz) in each session was 12 
and 3.3 min for PSR B0809+74 and B0950+08 
respectively. At low frequencies it was 35.5 min for 
PSR B0809+74 and 15.63 min for PSR B0950+08. 
Table 1 lists the time and frequency resolution 
{At and A/) and the used number of channels 
(iVch) at each frequency. Individual pulses of the 
pulsars were recorded on the computer disk in 
all channels with a period synchronized with the 
precomputed topocentric pulse arrival time. Then 
the signal in all channels was shifted in accordance 
with the dispersion shift at the given frequency 
(reduction to the highest frequency channel); in 
channels affected by interference (if any) the signal 
was replaced by an average value found from the 
adjacent channels. The amplification in all channels 
was reduced to the same value by normalization, 
so that the noise dispersion in all channels were 
equal to its average value found from channels not 
affected by interference. At frequencies 41-89 MHz 
the record was performed in a window 1.8 Pi (Pi 
is the pulsar period) and at 112 MHz in a window 
0.9 Pi. To improve the signal-to-noise ratio (S/N), 
we averaged individual pulses: from 15 pulses at 
112 MHz to 40 pulses at low frequencies. The 
averaging time T for both pulsars at all frequencies 
is listed in Table 1. 



3. CORRELATION ANALYSIS OF THE DATA 

To analyze the intensity variations of the pulse ra- 
diation as a function of frequency (channel number) 



within the bandwidth of the multichannel receiver, 
we formed spectra of individual pulses I{fk) by av- 
eraging the signal in a selected range of longitudes 
of the pulse and in a region outside the pulse, on 
the noise segment of the record (for the same num- 
ber of longitudes) iNiJk) for each channel. Here 
fk = f+{k — 1) A/ is the frequency of the kth chan- 
nel, and / is the frequency of observation. Then we 
subtracted the noise component lN{fk) from I{fk), 
and the spectra of individual pulses thus obtained 
were analyzed with the purpose to get information 
in the frequency as well as in the time domains. 
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Figure 1. (a) Mean pulse profile of PSR B0809-I-74 
at a frequency of 111.87 MHz as observed on April 
6, 2001. Vertical axis: intensity in relative (com- 
puter) units; horizontal axis: time inside the pulse 
in milliseconds. Longitude ranges (interval 1 and 
interval 2) in which the pulse intensity was aver- 
aged are marked, (b) Mean normalized CCF be- 
tween the spectra of adjacent pulses whose inten- 
sities were averaged in the above-mentioned longi- 
tude ranges: 1 (lower panel) and 2 (upper panel). 

Characteristic frequency (/dif) and temporal (tdif) 
scales of scintillations were determined using the cor- 
relation analysis. We computed the mean cross- 
correlation function (CCF) between the spectra of 
adjacent, noiseless pulses averaged in the given range 
of longitudes of the pulse emission; CCF was normal- 
ized to the product of rms deviations aia2. We con- 
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sidered a pulse as noiseless if its amplitude exceeded 
a level of 4(Tjv. The mean profile of the pulsar was 
derived by addition of all individual pulses in the 
given observational session. Since the pulse averag- 
ing time T (Table 1) was considerably shorter than 
the scintillation timescale, the decorrelation of the 
spectra of adjacent (averaged) pulses was insignifi- 
cant; in return, we completely eliminated the uncor- 
related component of noise. The frequency scale fd\f 
was determined as a frequency shift at which CCF 
decreased by a factor of two. 

Since we receive linearly polarized emission, we 
must take into consideration the influence of polar- 
ization on its frequency-time structure. At low fre- 
quencies the degree of polarization for pulsars stud- 
ied by us is high: about 60% for PSR B0809-h74 [13] 
at the longitudes of the leading part of the mean 
profile and much lower in its trailing part; (70- 
80)% for PSR B0950-f08 [14] throughout the pro- 
file. The rotation measure is RM — —11.7 rad/m^ 
for PSR B0809-f^74. The frequency of the Faraday 
rotation is 



PSR B0950+08 



Pf [MHz] = 17.475/^/RM, 



(1) 



where / is the observational frequency in hundreds 
of megahertz. Accordingly, Pp — 2140 kHz at 
112 MHz and Pp = 59 kHz at 41 MHz. The in- 
fiuence of polarization on the correlation analysis 
for PSR B0809-f74 is shown in Fig. 1. Longitude 
ranges in which we averaged the intensity to ob- 
tain the spectra used in the calculation of the av- 
erage CCF are shown. When averaging the lead- 
ing part of the pulse (interval i), a superposition 
of two structures— narrow-band (produced by prop- 
agation of the radiation in the interstellar plasma) 
and broadband (due to a much slower intensity vari- 
ation with frequency as a result of the rotation of the 
polarization plane across the receiver bandwidth)— 
is visible. When averaging the trailing part (inter- 
val 2), the slow component is virtually absent; this 
confirms the small degree of polarization at these 
longitudes. In the frequency-time data analysis of 
PSR B0809-f 74 we have used the average over the 
longitudes of the trailing part of the mean profile. 

For PSR B0950-I-08 the rotation measure is con- 
siderably smaller than for PSR B0809-I-74; however, 
in the literature its values appreciably differ: in the 
catalog of pulsars [15] RM = (1.35 ± 0.15) rad/m^, 
in [14] RM = (4-6) rad/m^, and in [16] RM = 
(-0.66 ± 0.04) rad/m^. Probably, authors [14] un- 
derestimated the contribution of the ionosphere to 
the obtained rotation measure. For RM = 1 rad/m^ 
the Faraday rotation frequency is Pp — 20.4 MHz at 
112 MHz and Pp = 1200 kHz at 41 MHz. The influ- 
ence of polarization on our observations at 41 MHz 
is shown in Fig. 2. The profile averaged over the 
session (Fig. 2a) is double-peaked with a separation 
between the components of 13 ms. Arrows shown 
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Figure 2. (a) Mean profile of PSR B0950-I-08 at a 
frequency of 41 MHz as observed on January 16, 2004. 
Vertical axis: intensity in relative (computer) units; 
horizontal axis: time inside the pulse in milliseconds. 
The time resolution is 2.56 ms, frequency resolution 
is 20 kHz. Arrows show longitudes at which the pulse 
intensity was selected in all channels for the calcu- 
lation of the correlation functions, (b) Mean CCF 
between the pulse spectra taken at longitudes 1-1, 
1-2, and 1-3. The curve in the lowermost graph 
is the fitted sine wave with a frequency of 600 kHz. 



longitudes at which we selected pulse intensities over 
all channels to form the spectrum at the given lon- 
gitude. Then we calculated mean normalized CCF 
between the spectra at longitudes 1-1 (autocorre- 
lation function, ACF), 1-2, and 1-3 for all noise- 
less pulses. Figure 2b shows mean CCF between the 
spectra at the corresponding longitudes. In these ob- 
servations we used the receiver with a channel band- 
width of 20 kHz. The narrow unresolved feature 
corresponding to the frequency scale of diffractive 
scintillations at 41 MHz has a maximum at the zero 
frequency shift, while the slow component is shifted 
to the left with an increase in the longitude offset of 
the spectra; this is due to the change in the position 
angle across the pulsar profile. Fitting a sine wave to 
the slow CCF component 1-3 yields the modulation 
frequency Pp = (600 ± 60) kHz, which corresponds 
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Figure 3. Mean normalized CCF for the 
spectra of adjacent pulses at 88.57 MHz (up- 
per graphs) and 62.43 MHz (lower graphs) for 
PSR B0809+74 (left) and PSR B0950+08 (right). 
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Figure 4. Mean profiles of PSR B0950+08 and 
PSR B0809-I-74 at 41 MHz as observed on Decem- 
ber 25, 2003 (left) and mean normalized CCF be- 
tween the spectra of adjacent pulses on this fre- 
quency (right). The time resolution is 5.12 ms, 
and the frequency resolution 1.25 kHz. In the 
left graphs, the vertical axis plots the intensity 
in relative (computer) units, and the horizontal 
axis the time inside the pulse in milliseconds. 



to RM = (2 ± 0.2) rad/m^. The observations were 
carried out in the nighttime, and we suppose that 
the contribution of the ionosphere to this value was 
not more than 1 rad/m^. Consequently, RM toward 
PSR B0950-h08 is about 1 rad/m^, and the effect 
of polarization on the frequency-time structure of 
the emission for this pulsar is considerably smaller 
than for PSR B0809-H74. Therefore, for obtaining 
the spectra we used the average over all longitudes 
determined by the level of 0.25 of the mean profile 
maximum of the pulsar. 

Figure 3 shows for both pulsars mean CCF be- 
tween the spectra of adjacent pulses at 88.57 MHz 
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Figure 5. Mean normalized CCF for the spec- 
tra of pulses taken at appropriate time inter- 
vals (horizontal axis); left: for PSR B0809-I-74 at 
111.87 MHz, right: for PSR B0950-I-08 at 88.57 MHz. 



Table 2. Characteristic scales of diffractive scintillations 



Parameter 


/, MHz 


41 


62.43 1 88.57 | 111.87 


/dif, kHz 

^dif, s 

/dif, kHz 
idif, s 


PS 
2 ±0.6 

PSI 

1.5 ±0.4 
350 ± 100 


R B0809+ 

7±3 

I B0950 + 

25 ± 10 
400 ± 100 


74 

20± 10 
450 ± 70 
08 

100 ± 40 
450 ± 120 


45 ±5 
600 ± 100 

220 ± 60 
>200 



(upper graphs) and 62.43 MHz (lower graphs) ob- 
tained in individual observational sessions. Table 2 
lists, together with their errors, the average values of 
characteristic scales of scintillations /dif found from 
all sessions. The errors correspond to standard de- 
viations from the mean. The slow CCF component 
is due to the effect of polarization. It is visible 
that with decreasing frequency the relative ampli- 
tude of the slow component increases; this testifies 
to an increase in the degree of polarization with 
frequency. Figure 4 presents the mean profiles of 
PSR B0809+74 and B0950+08 for one of the ob- 
servational sessions at the low frequency 41 MHz 
(left-hand graphs) and mean CCF between the spec- 
tra of adjacent pulses at this frequency (right-hand 
graphs). The frequency resolution here was 1.25 kHz 
for both pulsars. To improve the signal-to-noise ra- 
tio for PSR B0950+08 observed at / = 41 MHz with 
the narrow-band receiver, we used a time constant 
of 10 ms; therefore, the components of the mean 
profile merge. For PSR B0809+74 a well-resolved 
two-peak profile with a separation between compo- 
nents of 96 ms is observed. 

We determined the characteristic timescale of 
scintillations as the time shift at which the coeffi- 
cient of correlation between the spectra separated by 
an appropriate time interval decreases by a factor of 
two. Figure 5 presents the corresponding functions 
for both pulsars obtained from observations at one 
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of the frequencies: 111.87 MHz for PSR B0809+74 
(left) and 88.57 MHz for PSR B0950+08 (right). 
The average values of the characteristic scintilla- 
tion scales tdif found from all sessions are given to- 
gether with their errors in Table 2. The errors cor- 
respond to standard deviations from the mean. For 
PSR B0809-I-74 we did not manage to determine idif 
at 41 and 62 MHz because of a poor signal-to-noise 
ratio. 



4. STRUCTURAL ANALYSIS OF THE DATA 

We used the obtained frequency-time correlation 
functions for the structural analysis of the data and, 
accordingly, for the study of the spectrum of ISP in- 
homogeneities. The power-law spectrum of electron 
density fluctuations is deflned as 

$7Ve(9) = CL(Z-", (2) 

where characterizes the plasma turbulence 

level, q = ^n/p and p are, respectively, the spatial 
frequency and spatial scale of an inhomogeneity in 
the plane perpendicular to the Hne of sight (generally 
it is a three-dimensional vector) . The structure func- 
tion of phase fluctuations D{p) and spectrum 
are related by the Fourier transform. In the case of a 
statistically uniform distribution of inhomogeneities 
in the medium D{p) is described by relationship [17] 

D^ip) = (fceop)"-^ (3) 

(fcGo)"-' = A{n)iXrefClf,L/{n - 1), 

A( ^ ^ 24-"7r3 

[r2(n/2)sin(7rn/2)]' 

Here Qq is the characteristic scattering angle, A is 
the wavelength, k = 2tt/\ is the wavenumber, re 
is the electron classical radius, L is the effective dis- 
tance to the turbulent layer or distance to the pulsar 
in the case of a uniform distribution of turbulence 
along the line of sight. The spatial scale is related 
to the timescale by a simple relationship p = VAt, 
where V is the velocity of motion of the line of sight 
across the picture plane. If the pulsar velocity con- 
siderably exceeds the velocity of the observer and 
of the motion of the turbulent medium, then it de- 
termines the displacement of the line of sight. As 
was shown in [3], for small time shifts At the phase 
structure function Ds{At) can be obtained from the 
correlation function of intensity fluctuations Bi{At): 

D.(At) = M^I^?^ at At<U,, (4) 

where (J) is the mean intensity for an observational 
session. In the frequency domain we have used the 



relationship 

where A/ is the frequency shift. To reduce the data 
from different frequencies to a single frequency /o , it 
is necessary to scale temporal and frequency struc- 
ture functions in accordance with the law: 

£>,(/o, At(/o), A/(/o)) = D,{f,At,Af){f/fof. 

(6) 

As noted in [3], Af(/o) = Af(/); however, A/(/) ^ 
A/(/o). In the case of purely diffractive scintilla- 
tions we have 

A/d(/o) = ifo/ffAfif). (7) 

In the presence of a strong angular refraction, i.e., 
when the refraction angle considerably exceeds the 
scattering angle, the frequency dependence A/ is 
quite different: 

A/.(/o) = ifo/ffAfif). (8) 

We have used these formulas in the analysis of our 
data considering two models, diffractive and refrac- 
tive. 

Figures 6 and 7 show for both models frequency 
structure functions of phase fluctuations in the 
double logarithmic scale for PSR B0809+74 and 
PSR B0950+08 respectively. All the data were re- 
duced to the same frequency /o = 1000 MHz. As in 
our previous publication on the study of ISP spec- 
tra toward a number of pulsars [3-5], we have used 
this value of /o for convenience of a comprehensive 
analysis of all the data. In the flgures the data at 
different frequencies are shown by different symbols. 
In Fig. 6 a fllled circle marks a point obtained from 
observations at 408 MHz [18]. Statistical errors of 
the structure functions were estimated using equa- 
tion (B12) from [19]. From Fig. 6 it is visible that 
for PSR B0809+74 the diffractive model (Fig. 6b) 
describes the experimental data much better than 
the refractive model (Fig. 6a). A straight line least- 
square-fltted to the flrst points of the structure func- 
tions has a slope /3 = 1.3 ± 0.2. For fltting we must 
take the values at shifts smaller than or of the order 
of the characteristic frequency scale of scintillations. 
For PSR B0950-h08 (Fig. 7) the data at different fre- 
quencies in both models differ not so strongly as for 
PSR B0809+74; however, the scatter of the points 
for the refractive model (Fig. 7a) is smaller. The 
largest difference takes place for a point obtained 
from observations at 41 MHz (fllled triangle). The 
slope of the fltted straight line is /? = 0.96 ± 0.05. 

Figures 8 and 9 show temporal structure func- 
tions of phase fluctuations for PSR B0809-h74 and 
PSR B0950-I-08 respectively. In addition to our 
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Figure 6. Frequency structure functions of phase 
fluctuations for PSR B0809+74 in two models: re- 
fractive (a) and diffractive (b). Open circles: 111.87- 
MHz data; asterisks: 88.57 MHz; square: 62.43 MHz; 
triangle: 41 MHz; filled circle: 408 MHz [18]. The 
data for all frequencies have been reduced to the same 
frequency fo = 1000 MHz. The straight line least- 
square-fitted to the first points of structure functions 
for the diffractive model has a slope /3 = 1.3 ± 0.2. 
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Figure 7. Same as in Fig. 6, but for PSR B09.50+08. 
The fitted straight line has a slope 13 = 0.96 ± 0.05. 
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Figure 8. Temporal structure function of phase 
fiuctuations for PSR B0809-I-74. Open circles: 
111.87 MHz; asterisks: 88.57 MHz; small open 
squares: 933 MHz [19]. The data for all frequen- 
cies have been reduced to the same frequency fo = 
1000 MHz. The straight line has been least-square- 
fitted to the initial points of the structure func- 
tions. Top: axis of the corresponding spatial scales. 
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Figure 9. Temporal structure function of phase 
fluctuations for PSR B0950-I-08. Asterisks: 
88.57 MHz; small filled squares: 62.43 MHz. 
The data for all frequencies have been re- 
duced to the same frequency fo = 1000 MHz. 
The straight line has been least-square-fitted. 
Top: axis of the corresponding spatial scales. 



data, we have used in Fig. 8 also a structure func- 
tion obtained in [19] from observations at 933 MHz 
(open squares) . It is visible that the data at different 
frequencies for PSR B0809-I-74 presented in the dou- 
ble logarithmic scale are well described by a unified 
power-law spectrum. A straight fine fitted to the ex- 
perimental points has a slope a = 1.7 ± 0.04. As it 
is visible from Fig. 9, for PSR B0950-K08 temporal 
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structure functions at diflerent frequencies reduced 
to 1000 MHz satisfy a power law with an exponent 
a = 1.0 ± 0.05. The scale of the upper horizontal 
axis in Fig. 8 and 9 corresponds to spatial scales of 
interstellar plasma inhomogeneities p = V lS.t. Here 
we have used the pulsar velocities measured in [6]: 
V = 102 and 36.6 km/s for PSR B0809+74 and 
PSR B0950+08 respectively. 
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Figure 10. Characteristic frequency scale of scintil- 
lations as a function of the frequency of observation 
for PSR B0950+08 (left) and PSR B0809-I-74 (right). 
The filled squares show the scales obtained in this 
work, the asterisks the data of [20], the crust the data 
of [21], the open circle the data of ]22], the filled circle 
the data of ]10], the triangle the data of ]18], and the 
open squares: data from J7]. The lines are the result 
of least-squares fits. The slope for PSR B0950+08 us- 
ing our data is 7 = 4.9 ± 0.6, while the slope for the 
data at frequencies >112 MHz is 7 = 2 ± 0.6. The 
fitting for PSR B0809-I-74 fitting using all the points 
yielded the slope of 7 = 3.4 ± 0.15. The dashed line 
with a slope of 4.3 is the expected /dif(/) dependence. 
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Figure 11. Characteristic temporal scale of scin- 
tillations as a function of the frequency of observa- 
tion; top: PSR B0809-I-74, bottom: PSR B0950-I-08. 
Filled squares: scales obtained in this work, trian- 
gle: the data from J23], filled asterisks: J20], the 
open asterisk: ]24], the diamond the data of ]25], 
open circle: Jl9]. A straight line has been least- 
square-fitted to the data for PSR B0809-I-74; its 
slope is 1.1. Asymptotic straight lines with slopes 
of « = 2 and a = —0.5 corresponding to equa- 
tions (10) and (14) are shown for PSR B0950-f08. 



5. DISCUSSION 

As shown above, the temporal structure func- 
tions for PSR B0809-h74 and PSR B0950-H08 are 
described by a power law, and, accordingly, the 
spectrum of interstellar plasma inhomogeneities to- 
ward these pulsars is a power-law one with an ex- 
ponent n — a + 2; this corresponds io n = 3.0 ± 
0.05 for PSR B0950-h08 and n = 3.7 ± 0.04 for 
PSR B0809-I-74. Whereas for PSR B0809-I-74 the 
measured turbulence spectrum matches to the Kol- 
mogorov form (n = 3.67), for PSR B0950-h08 it 
is much flatter. In addition, we can say that to- 
ward PSR B0950-I-08 an appreciable angular refrac- 
tion takes place. As shown in [3], for the diffractive 
model the slope of the frequency structure function 
is related to that of the temporal structure function 
as (3 = a/2, and for the refractive model (3 = a. For 



PSR B0950-h08 we have (3 ^ a. For PSR B0809-f 74 
the accuracy of the determination of (3 is low, and 
within the 3ct error limits the relationship (3 = a/2\B 
fulfilled, though a more convincing evidence of the 
absence of angular refraction in this direction is a 
strong deviation of the data on the frequency struc- 
tural function from this model (Fig. 6). 

Figures 10 and 11 present the characteristic fre- 
quency and temporal scales of scintillations as func- 
tions of the frequency of observation for both pul- 
sars. Filled squares show our measurements, other 
symbols denote the data taken from [7, 10, 18- 
25]. For the reduction of the scales fd obtained 
by a different method in [21] to our measurements, 
we have used the relationship /jif — 0.3 fd- We 
have least-square-fitted straight fines to the obser- 
vation data. For PSR B0950-F08 (Fig. 10) fitting 
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was done separately to the low-frequency data of 
this work (slope 7 = 4.9 ± 0.6), and to the data 
at frequencies / > 112 MHz (slope 7 = 2 ± 0.6). 
For PSR B0809+74 fitting for all points (Fig. 10) 
yields 7 = 3.4 ± 0.15. As shown above, the fre- 
quency CCF of the flux fluctuations of the pulsar 
PSR B0809-I-74 is described by the model of diffrac- 
tive scintillations on inhomogeneities of the power- 
law spectrum of interstellar plasma turbulence with 
an exponent n = 3.7. For the diffractive model the 
characteristic frequency scale of scintillations is [17] 

/dif = c/(7ri?eo)2oc/2"/("-2), (9) 

where c is the velocity of light. For n = 3.7 the ex- 
pected for PSR B0809+74 value 7 = 4.3 is shown 
in Fig. 10 with a dotted line. This line is also con- 
sistent with the experimental points, except for the 
high-frequcncjf ones. It is possible that the values 
of /dif measured at high frequencies are underesti- 
mated because of the Hmited total frequency band 
of the observation. The characteristic timescale of 
the scintillations for the diffractive model is 

idif«iAeoFoc/2/("-2), (10) 

where k is the wavenumber and V is the veloc- 
ity of motion of the line of sight with respect to 
the turbulent medium. For n = 3.7 the expected 
for PSR B0809-h74 value 2/(n - 2) = 1.18; this 
agrees well with the experimentally obtained figure 
1.1 (Fig. 11). In [26] it was shown that the scin- 
tillation parameters for the pulsar PSR B0809+74 
and quasar B0917-I-624, which is located close to 
the pulsar on the celestial sphere, are determined 
by the same turbulent medium and this medium is 
distributed almost uniformly on a scale of the order 
of 500 pc. 

For the pulsar PSR B0950+08 the temporal and 
frequency correlation functions of the fiux fiuctua- 
tions are described by a more sophisticated model. 

Firstly, wc should take into account that the ob- 
servational data correspond to two modes of scintil- 
lations. At high frequencies (/ > scintillations 
arc weak, and at low frequencies (/ < fcr) the ob- 
servation data correspond to the diffractive compo- 
nent of strong scintillations. Near the critical fre- 
quency /cr scintillations correspond to the focusing 
mode, which has been poorly studied, either theo- 
retically or experimentally. According to our data, 
the scintillations index is about unity at all observa- 
tional frequencies; thus, scintillations are strong at 
frequencies / < 112 MHz. Consequently, the criti- 
cal frequency fa > 112 MHz. According to the data 
[25] obtained at / = 4.8 GHz, scintillations of the 
pulsar PSR B0950-I-08 are weak. In the same paper 
an estimate of the critical frequency < 700 MHz 
is given. In [27] on the basis of the compilation of 
the observational data of 1970s an estimate of the 
critical frequency ~ 300 MHz is given. 



Secondly, it is necessary to take into account the 
presence of a strong angular refraction. In this 
case, the characteristic frequency scale of diffractive 
(strong) scintillations is described by the relation- 
ship [3] 

/dif « c/TrrefOoeref CX /2+"/("-2)^ 
/ < /cr, 

where r^i is the effective distance to the turbulent 
layer, 8j-ef is the characteristic angle of refraction, 
60 is the characteristic scattering angle determined 
by (3). Forn = 3 equation (11) yields /jif oc f^. The 
observational data indeed correspond to a power law 
with a turnover. At low frequencies (/ < 112 MHz) 
the exponent 7 = 4.9 ± 0.6; this is consistent with 
the expected value for the mode of strong scintilla- 
tions in the presence of refraction. Therefore, we can 
assert that the critical frequency /cr > 112 MHz. In 
the case of weak scintillations we should replace the 
scattering angle ©o with the Fresnel angle 

Opr = {l/kr,ty/^. (12) 

As a result we obtain 

/dif W c/Trrefepreref OC /2-5. (13) 

At high frequencies / > 112 MHz the exponent 7 = 
2 ± 0.6; this also agrees with the expected value for 
the mode of weak scintillations. 

The characteristic timescale of scintillations for 
PSR B0950-I-08 is also described by a power law 
with a turnover. At low frequencies scintillations 
are strong, and the characteristic timescale of scin- 
tillations is determined by formula (10). At high 
frequencies scintillations are weak, and the charac- 
teristic timescale of scintillations is 

idif « (ref/A:)i/V^ef cx /"^Z^, (14) 

where V^f is the velocity of motion of the line of sight 
with respect to the turbulent layer. The asymp- 
totic relationships (10) and (14) are shown in Fig. 11 
with straight lines. We see that the theoretical de- 
pendences agree well enough with the experimental 
points. From the timescale of weak scintillations we 
can determine the distance to the layer responsible 
for the radiowavc modulation. If L is the distance 
from the observer to the turbulent layer and R is the 
distance from the observer to the pulsar, the effec- 
tive distance is 

1/ref = l/(ii - L) + 1/L. (15) 

The effective velocity is 

Ve( = Vo{R-L)/R + VpL/R, (16) 

where Vo and Vp are the components of the velocities 
of the observer and pulsar in the picture plane re- 
spectively. The pulsar velocity is = 36.6 km/s [6]; 
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therefore, Vp = Vq, and V^f = = 30 km/s. Using 
this value of Vet in (14) and supposing to = 15 min 
at the frequency / = 4.8 GHz [25], we obtain 

- L)/i? ^ 9.7 pc. (17) 

For this equation we find the solutions Li — 10 pc 
and Z/2 = 252 pc. Therefore, the turbulent layer is 
either near the observer or near the pulsar. It should 
be noted that parameters of the scintillations of the 
pulsar PSR B0950+08 arc rather close to those of 
the pulsar B0437— 47 [5], which is located at ap- 
proximately the same galactic longitude. For the 
latter pulsar it was shown in [5] that its scintillations 
as well as those of the quasar PKS 0405—385, at a 
small angular distance, are determined by a layer of 
the medium with enhanced turbulence at 10 pc from 
the observer. Therefore, we can accept that the pul- 
sar PSR B0950+08 scintillates on the same layer of 
the turbulent medium and that the right solution is 
Li = 10 pc. 



Table 3. Local turbulence levels 



PSR 


rad^ 


L, 
pc 


Ds{po)/L, 
rad^/pc 


AiVe, 

cm~^ 


B0329+54 


0.034 


1000 


3 X 10"" 


7.7 X 10-"* 


B0437-47 


0.0033 


10 


3 X lO^"* 


2.4 X 10"'^ 


B0809+74 


0.0003 


433 


1 X 10"'' 


1.4 X lO-* 


B0950-I-08 


0.002 


10 


2 X IQ-"* 


2 X 10-3 



Let us estimate the turbulence level toward 
PSR B0809-h74 and PSR B0950-h08. Note that 
we cannot use for the comparative analysis the val- 
ues of Cjfp, because they change their dimension 
with changing spectrum exponent n. We will use 
the values of the structure function on a selected 
scale po = 10^ m at a frequency / = 1 GHz. The 
local level of turbulence is related to the gradient 
Dsipo) as follows: {d/L)Ds{pa) « Ds{po)/L, where 
L is the characteristic thickness of the layer of the 
medium. The value of Ds{po)/L is proportional to 
the mean square of electron density fluctuations on 
inhomogeneities with a characteristic scale po; for 
Po = 10^ m we have: 

Ds{po)/L [radVpc] (18) 
- 50([AiVe(po)]') [ cm-6]. 

Table 3 Hsts the values of Ds{po) at a frequency 
/ = 1 GHz together with the values oi Ds{po)/L and 
ATVe for the pulsars studied by us in [3, 5] and in this 
paper. The distribution of the estimated local turbu- 
lence levels corresponds to three components of the 
interstellar medium mentioned in the Introduction. 
The greatest values of AA*"e correspond to the layer 
with an enhanced turbulence level (layer C): the pul- 
sars PSR BQ437-47 and PSR B0950-h08 scintillate 
on its inhomogeneities. The minimum value of AA^e 



obtained toward the pulsar PSR B0809-I-74 corre- 
sponds to a cavern. The intermediate value of AA^e 
is obtained toward the pulsar PSR B0329-I-54; it cor- 
responds to turbulent plasma in the space between 
spiral arms (component A). 



6. CONCLUSION 

As a result of the analysis of the phase struc- 
ture functions obtained by us from the study of in- 
terstellar scintillations toward PSR B0809-I-74 and 
PSR B0950-I-08 at low frequencies 41-112 MHz with 
invoking the previously published higher frequency 
data, we have shown that the spectrum of interstel- 
lar plasma inhomogeneities toward both pulsars is 
described by a power law. The spectrum exponent 
toward PSR B0950-I-08 differs appreciably from the 
Kolmogorov exponent; it is n = 3.00 ±0.05. Toward 
PSR B0809-(-74 the spectrum is a power law with an 
exponent n = 3.7 ±0.1; within the error limits this 
correspond to the Kolmogorov spectrum. The anal- 
ysis of the frequency dependence of the diffraction 
parameters has shown that toward PSR B0950+08 
the transition to the weak scintillation mode takes 
place approximately at a frequency of 200-300 MHz; 
here is a turnover in relationships idif(/) and /dif(/)- 
The experimental data match well enough the de- 
pendence expected by the theory. 

As a result of the analysis of the temporal 
and frequency phase structural functions for a 
purely diffractive model and model with a strong 
angular refraction, we have shown that toward 
PSR B0950±08 there exists a strong angular re- 
fraction. The conducted analysis of scintillations of 
this pulsar at different frequencies has shown that 
the distribution of inhomogeneities along the line 
of sight is not uniform and that the scintillations 
of PSR B0950+08 take place on a turbulent layer 
with enhanced electron density (layer C), which is 
located at a distance of about 10 pc from the ob- 
server. The same layer is responsible for the scin- 
tillations of PSR B0437-47 and and of the quasar 
PKS 0405—385. Here the local turbulence level on 
a scale po = lO"" m is 20 times higher than in the 
extended region with a size of the order of 430 pc 
responsible for the scintillations of PSR B0809+74. 
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